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■ Noise correlations, such as those observable in the time of flight images of a released cloud, are calculated for 
hard-core bosonic (HCB) atoms. We find that the standard mapping of HCB systems onto spin- 1/2 XY models 

, fails in application to computation of noise correlations due to the contribution of multiply occupied virtual 

Q . states in HCB systems. Such states do not exist in spin models. An interesting manifestation of such states 

(-7- ' is the breaking of particle-hole symmetry in HCB. We use noise correlations to explore quantum coherence of 

strongly correlated bosons in the fermionized regime with and without external parabolic confinement. Our 
0\ , analysis points to distinctive new experimental signatures of the Mott phase. 

(N ; 

In recent years, great experimental progress has been achieved in the coherent control of ultra-cold gases, which has been 
• used to attain laboratory demonstrations of iconic model quantum many-body systems. For example, by loading a Bose-Einstein 
condensate (BEC) into a tight two-dimensional optical lattice, an array of one-dimensional tubes has been created (JEHUS HI- 
Using this setup, recent experiments have been able to successfully enter the Tonks-Girardeau (TG) regime 0,01 predicted many 
• ' years ago @. Furthermore, by using an additional one-dimensional optical lattice in the direction of the tubes, a one-dimensional 
Mott insulator (MI) state with unit filling has been realized experimentally Q. In the TG limit, bosons behave like impenetrable 
particles - hard-core bosons (HCB) - and their dynamics are equivalent in many respects to those of a system of fermions. In this 
1 I paper, we calculate four-point correlation functions of HCB and the corresponding spin-1/2 and fermion systems. 

Experimental access to such correlation functions has recently been achieved 0|8|] by analysis of spatial (noise) correlations 
in shot-noise-limited images of expanded atomic clouds, following the original suggestion of Altman et al. @]. For example, 
for the case of bosons in a three-dimensional MI state, these noise correlations have been proven to complement the standard, 
first-order, characterization of phase coherence that is expressed by the direct image of the density of the expanded gas _ 
Such correlations have also been used to probe pair correlations of fermionic atoms generated by molecular dissociation | 

In this paper, we develop a theoretical framework to compute noise correlations for HCB systems confined in a one- 
dimensional lattice. The Hamiltonian of this HCB system is related to that of the spin-1/2 XY model, which in turn can be 
mapped onto the Hamiltonian of a spinless fermion system 

num. 

These three systems have identical spectra and local ob- 
servables, but their off-diagonal correlation functions differ. Experimentally relevant two-point correlation functions are identical 
j— 1 . for HCB and spin-1/2 systems, and the general formulation to calculate these correlations was developed by Lieb and Mattis 

■ mini, based on Wick's theorem. To our knowledge, explicit formulae to calculate higher order correlations have not been 
, worked out previously, except in certain asymptotic limits 1 19]. 

Here we calculate all four-point correlation functions and use them to compute noise correlations. One of the central results of 
this paper is the discovery of important differences between noise correlation functions of HCB and spin-1/2 systems. The root 
of this difference is the bosonic character of HCB that allows for multiple occupancy of the virtual states (MOV) even though 
the ground state is at most singly occupied. This is in contrast to spin-1/2 systems whose Hilbert space excludes multiple 
occupancy of all states. Here we show that the standard Lieb and Mattis lUOllllll formulation for spin 1/2 systems can still be 
utilized for calculating correlations functions in HCB systems provided the operators are written in normal order form by using 
, boson commutation rules, before applying any transformation. 

O ■ Our study shows that MOV states are responsible for various differences between the HCB and the XY spin-1/2 correlations. 
One particular way to characterize this difference is the particle-hole symmetry: in HCB noise correlations, this symmetry is 
broken while it is preserved in the corresponding spin systems. 

We use the noise correlations to characterize quantum coherence in both homogeneous lattice systems and systems subject to 
external parabolic confinement, as is often the case in experiments. Previous studies of noise correlations have been carried out 
in the MI limit our analysis extends this understanding into the strongly correlated (fermionized) regime. We show that, 

independent of the filling factor of the system, HCB exhibit second order coherence. It is manifested by peaks in the noise shot 
images which reflect the order induced by the lattice potential. This suggests that second order coherence in noise correlations is a 
generic attribute of correlated bosonic systems; as was also noted in earlier studies in a different context , such coherence 

does not merely signal reduced number fluctuations. Our detailed study involving all fillings factors highlights distinctions 
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between the Mott regime and the compressible regime. In the MI limit, the intensity of the density-density correlations between 
any two points separated by a reciprocal lattice vector is constant. In the compressible phase, on the other hand, the corresponding 
noise correlations depend strongly on the spatial position in the image. Thus a regular pattern in noise correlations could serve as 
a definite signature of the Mott phase. Additionally, the peaks of noise correlations in the compressible phase are accompanied 
by satellite dips immersed in a negative background. These dips are suppressed in the insulating phase. 

In Section I, we discuss various Hamiltonians related to HCBs and describe noise correlations that can be measured in exper- 
iments involving cold atoms in optical lattices. In Section II, we develop the theoretical framework needed to calculate all four 
point correlation functions relevant for computing noise correlations. In view of the relationships among HCB, XY spin-1/2, and 
spinless fermion models, we develop explicit formulas to calculate noise correlations for all of them. In Section III, we use the 
formulation of section II to study noise correlations for a spatially homogeneous gas of HCB. We discuss their basic properties 
and describe in detail how the strengths of the peaks and dips depend upon the filling factors. Our detailed numerical studies 
suggest that although the peak and dip amplitudes depend upon the size of the system, the peak to dip ratio is size independent. 
We compare noise correlations for HCB, spin-1/2 XY chains, and the corresponding fermionic systems, and show that the spin 
noise correlations preserve the particle-hole symmetry of the Hamiltonian, but that this symmetry is broken in HCBs. In Section 
IV, we compute the noise correlations for experimentally relevant parameters and compare results for Mott vs. compressible 
phases. Section V summarizes our results. 



I. INTRODUCTION AND FORMULATION 



A. Bose-Hubbard Hamiltonian and Observable 



The Bose-Hubbard Hamiltonian describes bosons in optical lattices when the lattice is loaded in such a wa y th at only the 
lowest vibrational level of each lattice site is occupied and tunneling occurs only between nearest-neighbor sites I13lf20ll : 

3 3 

Here aj is the bosonic annihilation operator of a particle at site j, hj = OjCij, and the sum (i, j) is over nearest neighbors. The 
hopping parameter J, and the on-site interaction energy U are functions of the lattice depth. Vj represents any other external 
potential such as a parabolic confinement or on-site disorder. 

In a typical experiment, atoms are released by turning off the external potentials at time t = 0. The atomic cloud expands, 
and is photographed after it enters the ballistic regime. Assuming that the atoms are noninteracting from the time of release, 
properties of the initial state can be inferred from the spatial images. As explained in detail in Ref. I9t ll4ll . the density distribution 
after the release, (h[x\) = (if!(x,t)^t(j(x,t)), with ip(x,t) the bosonic field operator, can be written in terms of the first band 
annihilation and destruction operators at lattice sites as J2i j w*(x — ia, t)w(x — ja, t){a\aj) where w(x,t) is the free evolution 
of a Wannier function centered around the origin and a is the lattice spacing. Assuming the Wannier functions can be described 
by Gaussians with initial width er /2, after time of flight t, the spatial density at position x, (n[x]), can be written as 



(n\x 



]) = w 2 (x, t) £ e^^'H"-™] (4 a m ) oc (n Q >, (2) 



|E^ M (4«J, (3) 



where w 2 (x, t) = y -^e 2x 1° and a = Td/(Ma ). The wave-vector P(x, t) = Mx/(ht) defines a correspondence 
between the position in the expanding image, x, and the lattice wave vector, Q. Here M is the particle mass and L is the number 
of lattice sites. 

Density-density correlations in the expanding cloud, Q[x, x'}, can also be linked to initial correlations in the lattice. For the 
evaluation of Q[x, x'] it is important to keep in mind that only normal ordered expectation values can be safely replaced by its 
projection into the lowest band. The second order correlation Q[x, x'] is given by H: 

G[x, x'} — 

(w(x)w(x')) 2 Y e ,p ( I *l"-™) e lP ( 1 ' t H ! -^tfl}fl J a m ) 

+w 2 (x)5(x - x'){h[x\) - (h[x]}(h[x']) (4) 
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where S(x — x') = LS x . x i is a delta function which arises from normal ordering in the continuum. Eq. @ can be mapped to 
quasi-momentum correlations in the lattice, Q[x, x'\ — ► A(Q, Q'), as: 

A(Q,Q') ee (h Q h QI ) - {h Q )(n Q >) 

~ (nQ)(5Q,Q>+nK - Sq,Q>), (5) 

where K = 2n/a is the reciprocal lattice vector and n is an integer. The first term in the second line of Eq. (jSJl arises due to the 
bosonic commutation relations of the operators and the identity 1/L YaLi e l2 ^ mnl / L = & m .nL- The quantity A(Q, Q') ee A qq i, 
that corresponds to density-density correlations in the expanding cloud of atoms, will be referred to as the noise correlations. 
The integer q characterizes the discrete values of the quasimomentum Q = j^q and this convention is going to be used hereafter. 

B. Hard Core Bosons (HCB) 
In the strongly correlated regime, Eq. (1) can be replaced by the HCB Hamiltonian lEoll . 

jj(HCB) = _ J^(gtg. +1 + St +i g.) + £ y. h . (6) 
3 3 

Here bj is the annihilation operator at the lattice site j which satisfies [bi^j, St] = 0, and the on-site conditions b 2 = bj = 0, 
which suppress multiple occupancy of lattice sites. The same relations are fulfilled by spin-1/2 raising and lowering operators, 
Oj~, and this is the reason why the HCB Hamiltonian is mapped to the spin 1/2 -XY Hamiltonian, 

= 2./ )►>:-;.: + a^ +1 ) + Vj £ ^ (7) 

3 3 

Here oi are Pauli matrices that satisfy anti-commutation relation at same sites. The parameter J describes the exchange inter- 
action between the XY-spins. A fully polarized state with all spins up corresponds to the Mott insulator state. We would like 
to emphasize that the mapping between the HCB and the spin Hamiltonian is a partial mapping as the two differ in the onsite 
commutation relation. 

For spin-1/2 systems, correlation functions are in general calculated by using the Jordan- Wigner transformation (JWT) fiol 
Eoll which maps spin operators into fermionic operators <t+ = cj exp — iri X)r =i cJjC„ . Here £j are fermion operators that 
obey anticommutation relations. This transformation maps the spin Hamiltonian into a non-interacting fermionic Hamiltonian, 

H(F) = - J E(sVi + + £ vjhj. (8) 

3 3 

Two-point HCB correlations have been calculated by invoking the HCB^ spin-1/2 — > fermion correspondence Jl7ll . As 
mentioned above, in contrast to the spin operators which obey fermionic commutation relations at the same site, HCB operators at 
equal sites obey bosonic commutation relations. Correlation functions such as (<fi\00' ! \4>) (where we have suppressed the lattice 
index for convenience ) are zero if O are spin lowering operators and non-zero if O are HCB field operators. The non-vanishing 
of these correlations can be traced to the fact that the HCB may have virtual states that have multiple occupancy: Recall that 
the HCB describes the U — ► 00 limit of the Bose-Hubbard Hamiltonian, and therefore, the limiting value of the Bose-Hubbard 
correlation functions must correspond to those of HCB. In general for bosons one can write \(f>) = A\0) + B\l) + e|2) + . . . 
where e ~ 1/U and A, B are constants. Therefore, ((f>\aa' ! \<fi) = A 2 + 2B 2 + 3e 2 + . . . . Taking the U — > 00 limit one gets 
(4>\aa' i \<fi) — > A 2 + 2B 2 > 0. On the other hand this correlation is always zero for the spin systems. 

The differences in the on-site commutation rules for bosons and spins have no effect on the calculations of local observables 
such as the density and the momentum distribution. However, correlation functions that involve a pair of operators bjVj, are 
affected by MOV states. In such cases, a direct application of the JWT is not correct. A simple recipe to take into account the 
MOV problem is to replace bjb^ by 1 + b^-bj before applying the JWT. We will refer to this recipe as the MOV rule. The validity 
of the MOV rule was checked by comparing various correlation functions obtained using the above recipe with those obtained 
by diagonalizing a full Bose-Hubbard Hamiltonian with a large U value (see below). 
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II. FOUR-POINT CORRELATIONS 



We now describe our calculations of four-point correlations. In view of the rather complex nature of the calculations, we will 
omit various technical details but focus on describing the final formulas and notations so that our results can be used for related 
future work. In addition to HCB, we will also discuss the related fermion and spin systems. 

For ideal fermions, four point correlation functions can be calculated using Wick's theorem directly. For example, the four- 
point correlations relevant for noise correlation are given by: 



(cJ,C m C z Cj) 



&lm9nj ~t~ QljQnm 9lm9nji 



where gi m are free-fermionic propagators: 



s=0 



(9) 



(10) 



with N being the total number of atoms and ipf' the s th eigenfunctions of the single-particle Hamiltonian — J{t/>j+i + 4-i) + 

Vi4 s) = e^4 s) . 

For HCB and spin systems, the calculation of four-point functions is more complex. Specifically, our calculations involve the 
following steps, (i) rearrange the operators so that the site index is ordered (this is only relevant for the case when three or more 
site indices are different); (ii) write operators in normal order form using bosonic rules to take into account the MOV; (iii) use 
the JWT accordingly to the prescription of Lieb and Mattis; (iv) use Wick's theorem to write higher order correlations in terms 
of two-point correlations. 

To present our results we denote the creation and the annihilation operators by b a , where a = +1(— 1) for annihilation 
(creation) operators, respectively, and label the "site ordered" four-point correlation function by Xabcd 



(11) 



In this equation it is implicit that a < b < c < d. This order is implied in all the expressions that follow. We define 
Gij = 2gij — 5ij and Bij = (b\ bj). The latter can be calculated in terms of GV,- lUlll . We introduce four matrices M, S, X 
and Y in terms of which our results for the correlation functions can be written as 



a/3 7 (5 _ Pi — 1 
Xabbd 2 

0/37(5 _ 1 — a/3 



X-abcd 



a/3-/S _ 1-7*5 



2 

T 



~(-l) sr 10 /l 

|M(a,M)| - f - + (1 - 8)^,-1 ) B ad 
\S(a,c,d)\+ (- + 
i(-l)"*|S( Cj a,&)| + (i + 



= (-1) 



b-\-d— c — a 



2 

- — "^ l X ( a ' b > c ' d )\ + 4 (^,-1 _ <5a,-i)l Y (a, b, c, d)\ 



(12) 
(13) 
(14) 
(15) 



M = 



X = 



/ G aa +i G a b-1 Gab+l ■■ G ac \ 

Gb-ia+i Gb-ib-i Gb-ib+i ■■ Gb-ic 

Gb+ia+i Gb+w-i Gb+ib+i ■■ Gb+ic 

\G C 

-la+l G c _ib_i G c -ib+l .. G c _ic / 

/ Gba ■■ Gbb-i Gbc+i ■■ Gbd-i Gbc \ 

G a +la ■■ Ga+lb-1 Ga+lc+1 ■■ G a +ld-l G a +lc 



Gb-la 
Gc+la 



Gb-ib-i Gb-ic+i 
Gc+ib-i G c +i c +i 



Gb-ld-l Gb-lc 
G c ,+ ld-l G c +i c 



\ G 



da 



G 



db-\ 



G 



dc+l 



G 



dd-l 



Gdc / 



(16) 



(17) 
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S = 



Y = 



G aa 
Gba 



Gab+1 
Gbb+1 



G a c \ 
Gbc 



V G c -i a G c -ib+l ■■ G c -i c / 

I G a -\-ia G a +lb Ga+la+1 ■■ G a +lb-l G a +lc+l 



Gb-la Gb-lb Gb-la+l 
G c +la G c +ib G c +i a +i 



Gd~la 
G C a 
V Gda 



Gd-lb Gd-la+1 
G c b G ca +i 
Gdb Gda+1 



Gb-ib-i 

Gc+lb-l 

Gd-lb-1 
Gcb-1 

Gdb-i 



Gb-ic+i 

Gc+lc+1 



Gd-ic+i 

Gcc+1 

Gdc+i 



(18) 



Ga+ld-l \ 

Gb-id-i 

Gc+ld-l 



Gd-ld-1 
Gcd-1 

Gdd-i 



(19) 



Here s is equal to for HCB and 1 for spins and rjp = <5g,-i- In the matrices dots indicate continuous variation of the indices. 
In the absence of dots, the indices are explicitly written and they may not change continuously. 

When three, four or two pairs of indices are equal the calculation of the four-point functions is simple and no determinant is 
actually required. The non-vanishing correlation functions of this type are given by: 



(blp n P n b n ) = g nn , (20) 

{bnb n bl n b rn ) = g n n9mm 9nmi (21) 

(P n b m b\ n b n ) = (-iy(9nngmm-9l m )+9nn, (22) 

(blbnblbm) = (b^bnblbn) = B nm . (23) 
with n 7^ m and Bij = (Upbj). The latter can be calculated in terms of Gy full . 



III. NOISE CORRELATIONS FOR HOMOGENEOUS SYSTEMS 



We now analyze various characteristics of the noise correlations that are calculated using the formulation discussed above. We 
consider particles in a periodic lattice in the absence of any external potential (Vj = 0). Here we compare and contrast HCB, XY 
spin-1/2 and the corresponding fermion system. This comparison illustrates the importance of MOV and the quantum statistics 
of the particles in noise correlations. Our study includes all filling factors v, < v < 1. Here v = N/L where L is the number 
of lattice sites and N is the total number of particles. We thus study the Mott regime (v = 1) as well as the non-Mott regime 
with 1 < v < 1. 

For the HCB and the spin systems the Toplitz-like determinants involved in the evaluation of noise correlations make the 
calculations complicated. The only simple case that can be studied analytically is when the system is a MI ( v = 1). In this case, 
9ji = an d me determinants involved in the calculation become trivial. The noise correlations are given by 

A HCB (Q u Q 2 )\,=i = S QltQ2 + \5 Q „ Q2+jK - |) , (24) 

where j is an integer. Eq. (I24> shows the characteristic bosonic peaks at every reciprocal lattice vectors, which reflect the order 
induced by the periodicity of the lattice. The peaks in HCB correlations reflect the bunching of the bosons due to Bose-Einstein 
statistics. It should be noted that the noise correlations at Qi = Q2 (autocorrelations) are different from those at Qi = Q2 + Kj, 
(j 0). Foraunit filled MI in particular, the former is two times larger than the latter ,A(<3i, Qi) = 2A(Qi,Qi + Kj). 

For spins, v = 1 corresponds to a fully polarized system. In this case, the local correlations dominate and the spins behave 
like fermions. The fermion correlations are given by 

A F (Qi,Q 2 ) = -n 2 Ql 5 Ql - Q2 , jK j^O, (25) 

where uq 1 = 1 for [Q\]k < Qf and zero otherwise. Here [}k means modulo reciprocal lattice vectors, j is an integer different 
from zero and Qp is the Fermi wave-vector Qp = 27r ^ V J 2 U (where || x || denotes the integer part). The fermionic noise 
interference pattern shows negative interference dips at every reciprocal lattice vector (except at Q\ = Q2 where A(Qi, Q\) = 
0) signaling the underlying anti-bunching due to Fermi-Dirac statistics. 
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From Eqs. 1241 and d25i one arrives at the conclusion that in the Mott regime {v = 1), HCB noise correlations with central 
peak height equal to 2 differ from the corresponding spin correlations which show no central peak. Thus this simple limit brings 
out the importance of MOV in HCB. 
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FIG. 1: Noise correlations at Q2 — for different values of Q — 2nq/(La) (In the plot the x-axis is in units of 1/a). The upper, middle 
and lower panels respectively correspond to HCB, spins and fermions (note the different scale used in the bottom panel). We use L = 55 and 

N = 27. 



For non-integer filling factors < v < 1, there is no simple analytic expression for the noise correlations. Therefore, all four 
point functions and the noise correlations are computed numerically. Our detailed investigation included lattices of various sizes 
and we studied both fixed end as well as periodic boundary conditions. Fig^illustrates the quantum coherence of HCB, spin 
and fermions as characterized by A q o- Our calculations show that for HCB as well as for spin systems noise correlations exhibit 
characteristic peaks at every reciprocal lattice vector for all filling factors. These peaks reflect the underlying order induced 
by the periodic lattice. The heights of the peak vary with the filling factor (see Fig. 2). The peaks are a manifestation of the 
boson nature of the atoms and are induced by interactions as they disappear in the non-interacting regime where all of the atoms 
are Bose-condensed. For fractional values of v, the peaks are accompanied by a small satellite dip at Q = 2ir/(La) + jK 
immersed in a negative background. In other words, the dips occur only in the compressible phase. These dips are present only 
for non-integer filling factor and are a manifestation of the quasi-long range coherence of the system. 

It turns out that Bogoliubov theory 11511 is helpful in gaining insight into the the dips and the negative background in the 
HCB noise correlations. The theory describes the weakly interacting regime and in principle can not be used to quantitatively 
describe the HCB gas. However, some aspects of the theory survive in HCB limit and suggest a possible mechanism for the dips 
and the negative background that characterize the compressible regime. The Bogoliubov approximation predicts the negative 
background A(0, Q) < and the satellite dip at low quasi-momenta (phonon regime). In the thermodynamic limit A(0, Q) 
diverges as — 4 ^ 2 . For finite lattices the height of the satellite dip does not diverge but scales like — ^j^j- Moreover in the 

Bogoliubov approximation the central peak scales like ^ cot [tt / L] so the peak to dip ratio remains finite (see also O). As 
we discuss below, the insensitivity of this ratio to the lattice size is also seen in HCB case. 

The negative background, A(0, Q) < 0, and the satellite dip at low quasi-momenta can also be qualitatively understood as 
follows: A(0, Q) can be written as A(0, Q) = (e|e) — NqUq where No is the number of condensate atoms, tiq is the number of 
atoms with quasi-momenta HQ (the population of which comprises the quantum depletion) and (e|e) is the amplitude of the state 
|e) = \a,Qdo\g), with \g) the many-body ground state. |e) represents the process of removing a particle from the condensate and 
one from the state with quasi-momentum HQ. In the absence of interactions all of the atoms are Bose condensed, correlations 
can be neglected and (e\e) = Nqtiq = 0. On the other hand, interactions dramatically change the behavior of the Bose gas. 
Collisions between zero-quasimomentum atoms admix into the condensate pairs of atoms at quasimomenta ±HQ. As a result 
pair excitations in the condensate become correlated so as to minimize the total energy of the system giving rise to the negative 
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sign in the noise correlations. For high quasi-momenta the interference plays a minor role but in the phonon regime n_Q diverges 
and the anti-correlations are maximal. 

As mentioned above for HCB bosons Bogoliubov theory is not qualitatively valid but it still gives a useful physical picture 
as the HCB ground state also exhibits off-diagonal quasi-long-range order that manifests in the power law decay of the density 
matrix and the macroscopic occupation of the zero quasi-momentum state (quasi-condensate) which scales like s/N fl7ll (see 
Fig. 5). 

In Fig. 2 we study the dependence of the peak and dip height with the filling factor for HCB and the correspondent spin 
system. As illustrated there the heights (Aoo) and the dips (| Aio|) of the peaks increase with the filling factor up to a maximum 
value and then begin to decrease. At v = 1, for HCB, the height approaches 2 while the dip vanishes (A(Q, 0) — > —2/L). 
However, for the spin model, the corresponding value is for both peak and dip. 




0.2 0.4 0.6 0.8 
V 



FIG. 2: The figure shows the central peak at occurs at q = q' = 0(top) and the the corresponding satellite dip intensity at q = 1 and q' = as 
we vary the filling factor. Red and blue respectively correspond to HCB and spin systems. Here L — 89 and we use free boundary conditions. 

The central peak attains its maximum value at v m 0.6 and the curve describing the variation of its height with the filling factor 
lacks reflection symmetry about v = 1/2. These results were confirmed for various lattice sizes as well as for periodic boundary 
conditions. The asymmetric behavior implies that the noise correlations for HCB do not preserve particle hole symmetry. In 
contrast, the satellite dip appears to preserve it. The asymmetry in HCB correlations has its roots in the MOV states as the 
spin correlations (which do not have MOV) preserves the particle-hole symmetry. In fact, the spin Hamiltonian, Eq. {7) with 
the additional constraints [ffZyjO^"] = 0, and [a~^ ,af]+ = 1, is particle-hole symmetric under the transformation a~ = h\, 
af = hi. The operators h\ and hi are the creation and annihilation operator for holes. On the other hand, if [<t~,<7 ? + ] + = 1 
is replaced by [hi, h\] = 1, as is the case for HCB, the particle-hole symmetry is no longer preserved under this transformation 
(which in terms of HCB operators reads as hi = h\, h\ = hi). 

It is important to mention that, in general, different observables do not have to preserve the symmetries of the Hamiltonian. For 
example, the momentum distribution for spins has an explicit dependence on the density and is not completely symmetric around 
v =1/2 lfl7ll . Nevertheless as shown in Fig. 2, noise-correlations do preserve the particle hole-symmetry of the Hamiltonian. 

Our detailed numerics shows that the heights of the peaks and the dips depend upon the lattice size and are expected to diverge 
in the thermodynamic limit. However, our simulations also suggest that the ratio of height to dip is size independent (behavior 
also observed in the Bogoliubov calculations) for various filling factors. In particular for HCB the ratio near the maxima of the 
central peak is found to approach 2tt (See Fig. 3). The scale invariance of this ratio makes it a more desirable physical quantity 
that can be compared in different experiments. 

The striking differences between the spins and the HCB noise correlations not only highlight the importance of including 
MOV states but also explicitly show that a direct mapping between the two systems is not valid for all observables. We checked 
the validity of MOV rule by calculating noise correlations using a direct diagonalization of the Bose-Hubbard Hamiltonian for 
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small size systems. Fig. 4 illustrates the variation of height of the central peak and the dip as a function of the interaction U. As 
U increases noise correlations reach the asymptotic value that was found to be in excellent agreement with the value calculated 
from the HCB model including MOV states. We also calculate the noise correlations for the spin model, to emphasize the fact 
that MOV corrections can be significant. Even for a small size system, Fig. 4 clearly shows broken symmetry for HCB while the 
symmetry is preserved for the spin system. In Fig. 4 the central peak is found to be associated with a dip for all values of U in 
consistency with the Bogoliubov calculations. 
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FIG. 3: The upper and the lower panels respectively show the ratio of height to dip for the HCB and spin systems. The three curves in each 
panel respectively correspond to L = 89 (dotted), L = 55 (short dashed), L = 21 (long dashed) lines. 



IV. MOTT VS NON-MOTT PHASE IN A PARABOLICALLY CONFINED CASE 



In this section we study noise correlations for the experimentally relevant case when there is an additional parabolic confine- 
ment Vj = Vlj 2 . The advantage of having the parabolic confinement is that in this case it is always possible to realize a unit MI 
at the trap center for any number of particles with an appropriate choice of the trapping potentials I17lll8ll . In the homogeneous 
case on the contrary only for the integer filled case is a Mlpossible. In our analysis, the chosen parameters correspond to typical 
experimental set-ups such as the ones reported in Ref. Q|, We present our results for two different values of 51/ J, 0.008 and 
0.17. 

The density profiles (Fig. 5 top panel) show that for the case il/J = 0.17 the ground state of the system corresponds to a 
MI with all the central N sites singly occupied. For ft/ J = 0.008, on the other hand all, the sites have filling factor less than 
unity. The formation of a MI state with reduced number fluctuations and the localization of the atoms at individual lattice sites 
are signaled by the momentum distribution (Fig. 5 bottom panel), which shows a flat profile for ft/ J = 0.17. On the other hand, 
the non-vanishing off diagonal coherence and large number fluctuations in the fi/J = 0.008 case are reflected by the sharp peak 
in the momentum distribution at Q = 0. 

In Fig. 6 we plot A(Qi,Q2) for the Mott and non-Mott systems as a function of Qi and Q2. Noise correlations show 
interference peaks when Q\ = Q2 + nK. The existence of these peaks in MI as well as in non-Mott phase implies that the 
presence of peaks is not a signature of the insulating phase. Our numerical calculations show that, nevertheless, information 
about the insulator or superfluid character of the system can be extracted from noise-interferometry as only when the system is 
a MI, the noise-correlations exhibit a regular pattern, i.e. A(Qi, Q2) ~ — Q2) (small differences seen in Fig. 6 are due 

to the finite trap). In contrast, when tunneling is allowed and non-local correlations are established, they cause a decay in the 
intensity of the peaks as one moves away from [Qi]k = [Q2] k = 0. This decay is similar to the one observed in the momentum 
distribution away from Q = and reflects the non-vanishing off-diagonal coherence of the system. 
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FIG. 4: The figure shows the peak (upper) and the dip(lower) intensities for L — 6 and N equal to 2(dashed blue line) and 4 (dotted- 
dashed black) as the interaction U varies. The solid horizontal red lines show the corresponding result for HCB obtained using the theoretical 
formulation described in section I. The green( thick) line shows the result without MOV contribution and hence describes the noise correlations 
for the spin- 1/2 XY chain. In this case, v — 1/3 and 2/3 results coincide reflecting particle-hole symmetry. 



0.8 
0.6 




-20 -10 10 20 
J 

4 i 




' ■ ■ ■ ■ ■ ■ 1 

n -tt/2 n/2 n 

Q 



FIG. 5: Top: Density, Bottom: momentum distribution (Q is in units of 1/a) for different trapping potentials: 0.008(blue solid line ), and 0.17 
(red dashed-line)). Here TV = 19 and L — 55. The correlation functions are renormalized by a scaling factor N/Z where Z are the number of 
sites with non-zero density. 
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FIG. 6: Noise correlations as a function of Qi and Q2 for the same parameters used in Fig.4. Qi and Q2 are in units of 1/a. The upper and 
lower correspond to fl/J = 0.008( Non MI) and 0.17(MI) respectively. The correlation functions are renormalized by a scaling factor N/Z 
where Z are the number of sites with non-zero density. 



Furthermore, Fig. 6 shows that the satellite dips observed in partially filled translationally invariant systems are still present in 
the inhomogeneous case. When 0/ J = 0.008 they are located at = 2ir/(La) [Q2\k = and disappear in the MI phase. 

All of these results suggest that noise interferometry provides valuable complementary information about the phase coherence 
of the system and thus can be a suitable experimental tool for characterizing the MI phase. Experimental techniques capable of 
thoroughly characterizing the MI phase are crucial for the neutral atom based quantum information proposals. 

V. SUMMARY 

In this paper, we have derived explicit formulae to calculate four-point functions in HCB, spin-1/2 and fermionic systems 
in the presence of an external potential. Our results are valid for all filling factors and hence include the Mott as well as the 
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non-Mott regimes. Our work generalizes the Lieb and Mattis formulation for two point correlations to four-point correlations 
and is therefore an asset for all future studies of higher order correlations in HCB as well as in spin-1/2 systems. 

One of the key results of our studies is the fact that although multiply occupied states are suppressed in the ground state 
of strongly correlated bosons, MOV states have to be included for a proper calculation of correlations. These states lead to 
differences between HCB and spin-1/2 XY chains which have important manifestations such as the breaking of particle-hole 
symmetry in the HCB systems. Recently, there have been various proposals to use bosonic atoms in optical lattices as quantum 
simulators of spins models l2lll . Our study points out that such mappings have to be done very carefully as MOV states can 
certainly change the spin character of certain observables. 

We also showed that the noise correlation pattern is a sensitive probe to characterize the MI phase as: i) only in the MI phase 
is the noise-pattern completely regular, ii) there are additional satellite dips accompanying the Bragg interference peaks in the 
absence of a MI. These findings might have some impact in current efforts to implement a quantum computer using optical 
lattices as most of proposals require a well characterized MI state to initialize the quantum register filjll . 

We would like to emphasize that we have studied experimentally measurable shot noise correlations for strongly correlated 
bosons in the TG regime. Discussions related to the corresponding spin system is included to highlight the importance of 
multiple occupancy of virtual states in HCB, as the correlations in the spin systems describe correlations of HCB without MOV. 
In other words, results for spins are not presented in the context of any experiments in these systems. However, our general 
expressions for four-point correlations for spins will potentially be very useful for various studies of higher order correlations in 
spin systems such as XY and Ising chains with or without disorder. 
Acknowledgments 

We are grateful to Trey Porto for various suggestions and comments. We would like to thank L. Mathey et. al fl^l for sharing 
their results (under preparation) of a related study of noise correlations using bosonization techniques where a characteristic 
peak, followed by a satellite dip (similar to those shown in Fig.l) are found. This work is supported in part by the Advanced 
Research and Development Activity (ARDA) contract and the U.S. National Science Foundation through a grant PHY-0 100767. 
A.M.R. acknowledges additional support by a grant from the Institute of Theoretical, Atomic, Molecular and Optical Physics at 
Harvard University and Smithsonian Astrophysical observatory. 



[1] Laburthe Tolra B et al 2004 Phys. Rev. Lett. 92 190401 . 

[2] Moritz H, Stoferle T, Kohl M and Esslinger T 2003 Phys. Rev. Lett. 91 250402. 
[3] Paredes B et al 2004 Nature 429 277. 

[4] Kinoshita T, Wenger T and Weiss D S 2004 Science 305 1125. 

[5] Fertig C D et al 2005 Phys. Rev Lett. 94 120403. 

[6] Girardeau M 1960 J. Math. Phys. 1 516 . 

[7] Greiner M et al 2005 Phys. Rev. Lett. 94 110401. 

[8] Foelling S et al 2005 Nature 434 481. 

[9] Altman E 2004 et al, Phys. Rev. A 70 013603. 
[10] Lieb E H and Liniger W 1963 Phys. Rev. 130 1605 . 
[11] Lieb E , Schultz T and Mattis D 1961 Ann. Phys. 16 407. 
[12] Mathey L, Altman E and Vishwanath A 2005 cond-mat 70507 108 . 
[13] Jaksch D et al 1998 Phys. Rev. Lett.81 3108. 
[14] Roth R and Burnett K 2003 Phys. Rev. A 67031602. 
[15] Bogoliubov N N 1947 J. Phys 11, 23. 

[16] Similar results using bosonization techniques for ID systems have been found by Mathey L, Altman E and Vishwanath A . In preparation. 
[17] Rigol M and Muramatsu M 2005 Phys. Rev. A 72 013604. 

[18] Pupillo G, Rey AM, Brennen G B , Williams C J , and Clark C W 2004 J. Mod. Opt. 51 2395; Brennen G B , Pupillo G, Rey A M, Clark 

C W and Williams CJ 2005 J. Phys. B 38 1687. 
[19] Korepin V E , Bogoliubov N M , and Izergin A G 1 993 Quantum Inverse Scattering Methods and Correlation Functions, ( Cambridge 

University Press). 

[20] Sachdev S 1999 Quantum Phase Transitions ( Cambridge University Press). 
[21 ] Duan L M , Dernier E and Lukin M D 2003 Phys Rev Lett. 91 090402. 



